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A theoretical and experimental investigation was made of a porous 
wall diffuser used with a low density hypersonic nozzle. The Reynolds 
number range of the experiment varied from 1000 to 20,000 based on the 
nozzle diameter. At the low Reynolds numbers, nearly all of the flow 
passed through the pores of the diffuser. At the higher Reynolds num- 
bers, 70 to 85 percent of the flow passed through the throat of the 
diffuser. The measured pressure recoveries varied from 1 to 10 times 
the test section normal shock pressure. When models were introduced 
into the test section stream, the mass flow and pressure recovery of 
the diffuser were markedly reduced. Although the model used to describe 
the flow through the porous wall appears incorrect, the theoretical and 
experimental pressure recoveries and mass flows were in good agreement. 
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TECHNICAL MEHORANDIM X- 5 3410 


STUDY OF POROUS WALL LOW DENSITY WIND TUNNEL DIFFUSERS 

SUMMARY 


A theoretical and experimental investigation was made of a porous 
wall diffuser used with a low density hypersonic nozzle. The Reynolds 
number range of the experiment varied from 1000 to 20,000 based on the 
nozzle diameter. At the low Reynolds numbers, nearly all of the flow 
passed through the pores of the diffuser. At the higher Reynolds num- 
bers, 70 to 85 percent of the flow passed through the throat of the dif- 
fuser. The measured pressure recoveries varied from 1 to 10 times the 
test section normal shock pressure. When models were introduced into 
the test section stream, the mass flow and pressure recovery of the 
diffuser were markedly reduced. Although the model used to describe the 
flow through the porous wall appears incorrect, the theoretical and 
experimental pressure recoveries and mass flows were in good agreement. 


I. INTRODUCTION 


The near-free molecular flow regime represents a low density flow 
regime of great practical interest. In this regime, the flow about a 
body is determined by the character of both intermolecular collisions 
and molecule-wall collisions. The relative importance of these two 
types of collisions is primarily dependent upon the Knudsen number (Kn) 
which is the ratio of the mean free path to the pertinent test object 
dimension. As the Knudsen number increases, a smaller fraction of 
molecules that rebound from test object will collide with molecules 
that are proceeding to the test object, so the incoming scream of mole- 
cules becomes less and less dependent upon the test object. In the 
limit Kn the incoming stream is independent of the test object, and 

the flow is completely free molecular. 

To investigate this flow regime experimentally, it is necessary to 
provide a facility capable of producing a range of mean free paths that 
are greater than and less than the test object dimensions. Since the 
typical low density hypersonic boundary layer is one or two orders of 
magnitude greater than the mean free path, it is necessary that the wind 
tunnel be two or three orders of magnitude larger than the test object to 
obtain reliable results at high Knudsen numbers. Because of this large 
ratio, the wind tunnel must be large even for models with dimensions as 
small as one inch. 



This requirement for a large tunnel diameter means that the tunnel 
pumping system must also be large in terms of volume flow. This results 
from the large physical size of the tunnel and the presence of the low 
momentum flow in the boundary layer. In a typical low density diffuser, 
the relatively thick boundary layer prevents the attainment of any signi- 
ficant pressure recovery in terms of the test section normal shock pres- 
sure recovery through a diffuser. 

Bottorff and Rogers (1963) showed that a nozzle with porous walls 
could be used to provide some control of the boundary layer height through 
boundary layer suction. It was found that the boundary layer suction 
reduced the thickness of the boundary layer and therefore allowed a reduc- 
tion of the physical size of the nozzle for specified test section condi- 
tions . 

Previously preliminary experiments by Rogers (1962) showed that a 
porous nozzle and a porous diffuser could be used to obtain pressure 
recovery in excess of test section normal shock pressure recovery. Since 
these levels of pressure recovery would permit one or two orders of magni- 
tude reduction in the pumping speed requirement for a low density wind 
tunnel, it appeared that this approach might make it economically feasible 
to develop large low density facilities. The present report is a study 
of the operational characteristics of porous wall low density diffusers 
when operated with a low density wind tunnel in the Reynolds number range 
of 10^ to 1.2 X 10^. 


II. THEORETICAL CONSIDERATIONS OF DIFFUSER PERFORMANCE 


A. Nonporous Nozzle and Diffuser Walls 

In the typical operation of a supersonic wind tunnel (Figure 1), 
air is expanded from the nozzle throat (Station *) to the test section 
(Station 1) and compressed from the test section to the diffuser throat 
(Station 2). The flow becomes subsonic at Station 3. If the flow 
between the nozzle throat and the diffuser throat were completely isen- 
tropic, it would be possible to compress the flow to the sonic condition 
at the diffuser throat. Assuming the flow was adiabatic, the diffuser 
throat would be the same size as the nozzle throat. There would be no 
loss in stagnation pressure through the wind tunnel, and this condition 
would have the minimum possible pumping requirements. In practice, the 
flow is never isentropic between the nozzle throat and the diffuser throat, 
and frequently it is not adiabatic. Boundary layer flows and shock waves 
are the two main phenomena that prevent the flow from being isentropic. 
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The decrease in stagnation pressure resulting from the nonisen- 
tropic flow requires that the diffuser throat be larger than the nozzle 
throat. From the other viewpoint, the maximum attainable area ratio 
between the nozzle throat and the diffuser throat can be used as a 
measure of the stagnation pressure loss between the two throats. Assum- 
ing that the flow at each throat is uniform and one-dimensional, the con- 
tinuity equation can be used to develop the relationship between the 
throat conditions. The present development is simplified by considering 
the case of a nozzle and a diffuser having solid walls so that mass flow 
entering the nozzle throat must pass through the diffuser throat: 


P2 U2 As = % A^, 


( 1 ) 


but 


pu = — M = 


tI 


y + 1 

/ y . 1 a"2(7 -1) _ 

f 1 + ^-2 ■ j M 


( 2 ) 


Therefore, neglecting possible changes in 7 and R, 


7~H 


-. 2 ( 7 - 1 ) 



The tern (A/Aj.) is the area ratio associated with M2. 

tC 2 

The maximum attainable value of A-k/k^ corresponds to choking, so 
that M2 = 1 and (A/A*)^ = 1 . In this limiting case, the pressure recovery 
is given by 
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( 5 ) 



= — 4T7t 

A ^ •‘■rn t 


If the flow is adiabatic, pressure recovery at a choked 

diffuser throat is directly proportional to the area ratio between the 
nozzle throat and the diffuser throat. 

The result corresponds to the minimum pressure recovery that can 
be obtained at the specified area ratio. This follows from the fact that 
the sonic condition corresponds to the maximum mass flow per unit area. 
Therefore, with a fixed area ratio, a reduction in pressure recovery would 
have to be accompanied by a reduction in mass flow. Since the mass flow 
is fixed by the continuity equation, it is not possible to have a pressure 
recovery at the throat less than that given by equation (5). There may be 
significant losses downstream of the throat, so that the overall pressure 
recovery (Px^/^^T ) below that given in equation (5). This is illus- 

trated in Figure 2, which compares the results of equation (5) with the 
experimental results reported by Johnston and Witcofski (I960). The 
experimental pressure recovery exceeds the theoretical value for smaller 
values of area ratio but as the maximum experimental values of 

A^/As are approached, the experimental pressure recoveries fall below the 
theoretical ones. This is the result of losses downstream of the diffuser 
throat section. 

Figure 2 also emphasizes that experimentally the point of maxi- 
mum pressure recovery does not correspond to a choked throat. This 
follows from the prior argument that the maximum mass flow per unit area 
corresponds to the choking point. Contraction beyond this choking point 
combined with decreasing pressure recovery is not possible without corre- 
sponding decrease in mass flow. Since the mass flow was cpnstant during 
the experiments, it follows that the throat was not sonic at the maximum 
pressure recovery point. 

While the analysis based on a sonic diffuser throat is not valid 
in calculating the overall pressure recovery, it is useful in establishing 
the minimum pressure recovery that can exist at the diffuser throat. It 
is also possible to use the one-dimensional analysis to establish the 
maximum pressure recovery that can be obtained for a specified area ratio 
between the nozzle throat and the diffuser throat. This is done by assum- 
ing that the flow is isentropic between the two throats. In this case, 
the flow enters the diffuser throat at a supersonic velocity, and the 
idealized pressure recovery will be that associated with a normal shock 
at this supersonic diffuser throat Mach number. It can be seen that this 
is the maximum possible pressure recovery for the specified area ratio by 
observing that, while any consideration of boundary layer displacement 
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effects tends to reduce the Mach number, and therefore the shock losses 
associated with the flow, the viscous losses associated with the boundary 
layer growth are greater than the benefit resulting from the lower Mach 
number . 


Equation (4) can be used to calculate the maximum possible pres- 
sure recovery. In the equation, M 2 is the Mach number downstream of the 
normal shock wave. In the limiting case of a very high Mach number 
upstream of the normal shock, the downstream Mach number becomes 

^^upstream 

When equation (6) is substituted into equation (4) and the results evalu- 
ated for 7 = 1.4, 




/Tt . 

2 


( 7 ) 


While this equation is based on the assumption of a very high incoming 
Mach number, it is within 10 percent for throat Mach numbers down to 
M = 4.2. Now equations (5) and (7) represent two limiting sets of assump- 
tions. Equation (5) is based on the assumption that the stagnation pres- 
sure losses are such that the diffuser throat is choked and therefore 
represents the minimum pressure recovery that can occur at that throat 
area ratio and temperature ratio. Since equation (7) is based on the 
assumption that the stagnation pressure is constant to the diffuser throat 
where the flow becomes subsonic through a normal shock, it represents the 
uiaximum pressure recovery that can occur at that throat area and tempera- 
ture ratio. Since equations (5) and (7) differ only by the factor 1.65, 
it is clear that the pressure recovery is essentially established by the 
throat area ratio and the temperature ratio, and is independent of the 
processes between the two throats. The processes between the two throats 
will of course determine the possible throat area ratio and temperature 
ratio, but the significant point is that if the limiting values of throat 
area ratio and temperature ratio are established, the maximum throat pres- 
sure recovery must fall within the rang'e of equations (5) and (7). It is 
to be emphasized that the minimum measured pressure recovery may be below 
the limits established by equation (5). This is due to losses occurring 
downstream of the throat. 
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B, Porous Wall Nozzle and Diffuser 


It is necessary to determine the effect of the flow through the 
porous walls when calculating the pressure recovery of a system having 
porous walls. Since both the nozzle and diffuser walls may be porous, 
the decrease in mass flow between the nozzle throat and the diffuser 
throat must be considered in deriving the porous wall counterparts of 
equations (5) and (7). 

If the mass flow through the diffuser throat is less than the 
mass flow through the nozzle throat, equation (1) can be rewritten 


A - a 

P2 U2 As - ^ p* u* A^. 

•k 

This is the only change due to the suction. The porous wall counterparts 
of equations (5) and (7) become equations (9) and (10), respectively. 




* min 




(9) 


^ T*-^max 


A* As Ts 


( 10 ) 


Thus, equations (9) and (10) are just equations (5) and (7) multiplied 
by ifis/m^. 


To use these equations in calculating the limits on pressure 
recovery as a function of area ratio, it is necessary to determine the 
variation of m^/m^ and Tx^/^t.^ with area ratio. 

The term m 2 /m^ can be obtained if the suction flow through the 
porous wall can be determined. 



m 


suction 

m 

* 


( 11 ) 
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The temperature ratio T'p /T.^^ cannot be determined in such a 
direct manner. However, since the pressure recovery varies as the 
temperature ratio to the one-half power, the accuracy of the temperature 
ratio is less significant than the accuracy of the mass flow ratio. The 
temperature ratio can be estimated from the heat transfer to the walls of 
the nozzle and the diffuser by the following heat balance. 


. r- 


C = / h(T ^ - - T dA + (m^ - m 2 ) C T 

^ P J ad-wall wall^ * P suction flow 

A 


Heat In Heat transferred to wall 


Heat carried through 
porous wall 


( 12 ) 


”2 Cp T^^ 

Heat passing through 
diffuser throat 


This will be simplified by the following assumptions: 

(1) The Prandtl number is unity so that Tt = T 

(2) the wall temperature is constant = T^, and 

(3) the gas leaves the porous wall at the wall temperature 

so that T = T T T = T . 

suction wall w 

Using these assumptions, the temperature ratio becomes 



To solve this, it is necessary to determine the heat transfer coefficient 

(h). 
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The Reynolds analogy = KCj can be applied to the present 
analysis. If the effect of pressure gradient is neglected, the con- 
stant K is equal to one-half for a suction boundary layer as well as 
for a nonsuction boundary layer. As will be seen later, the effect of 
pressure gradient is small as far as the flow properties are concerned; 
thus, this effect will be ignored in calculating the heat transfer rate. 

The variation of pressure recovery with area ratio can now be 
estimated if the local skin friction coefficient and the total mass flow 
through the pores can be determined. Bottorff and Rogers (1963) pre- 
sented a technique for calculating the flow properties in a porous 
laminar flow nozzle. This technique used an iterative procedure which 
satisfied the momentum and continuity equation at each station in the 
nozzle. The energy equation was incorporated by assuming a Prandtl 
number of unity and using the Crocco integral relationship to relate 
momentum and energy (Crocco, 1948). The local boundary layer charac- 
teristics were determined from the suction boundary layer calculations 
of Iglisch (1949) modified to include the effects of heat transfer, com- 
pressibility, and axis 3rtnrae trie boundaries. These characteris tics were 
used to calculate the wall shearing stress, the momentum thickness, the 
displacement thickness, and the height of the boundary layer. 

The flow through the pores was estimated by assuming that the 
pores were connected to a reservoir containing a gas at uniform condi- 
tions. The stagnation pressure for the reservoir gas was assumed equal 
to the local static pressure in the nozzle, and the stagnation tempera- 
ture of the reservoir gas was assumed equal to the nozzle wall tempera- 
ture. The effect of pore L/D and Reynolds number was included in the 
calculation of the pore mass flow. Since this approach provides the 
flow characteristics necessary to calculate ma/^* 

nozzle, the computer program was modified to include the calculation of 
the flow into the diffuser. 


C. Starting Limitations 

Since the pressure recovery is directly proportional to the 
ratio of the nozzle throat area to the diffuser throat area, it is 
important to be able to establish a large value of this ratio. In a 
solid wall diffuser, the value of this ratio is limited by the starting 
process. During the starting of a supersonic nozzle, the shock waves 
must pass through the test section into the diffuser throat, which must 
be sized to allow all of the relatively low pressure recovery air to 
pass through the diffuser throat during the starting process. This 
places a severe limitation on the pressure recovery of a fixed geometry 
diffuser, and to obtain larger pressure recoveries, it is necessary to 
use variable geometry diffusers. Since this is not practical for an 
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axisymmetric diffuser, most axisymmetric diffusers have fixed geometry 
and limited contraction* While the porous diffuser must operate under 
a similar limitation, it is in a more favorable position. This results 
from the fact that the pores, as well as the diffuser throat, are avail- 
able for the passage of mass flow during the starting process* Thus, 
during the critical phase when the starting shock system is near the test 
section, the entire porous diffuser and throat are available for mass 
flow, whereas the solid wall diffuser has only the throat area available. 
After the porous wall diffuser has started, the low static pressure com- 
bined with the supersonic flow field acts to reduce the flow through 
the pores so that, whereas the entire flow could pass through the pores 
with the shock system at the test section, after the diffuser has started, 
only a fraction of the flow passes through the pores. 

To calculate the mass flow that passes through the pores of the 
diffuser during the starting process, it is necessary to make several 
assumptions* These will be discussed in the section on the results of 
the numerical calculations . 


III. THEORETICAL CALCULATION OF DIFFUSER BOUNDARY LAYERS. 


A* General Approach 

The method of Bottorff and Rogers (1963) for calculation of 
compressible laminar boundary layers in nozzles with suction was extended 
during this study to include calculation of the flow in porous diffusers. 
The extension assumes that the boundary layer concept can be carried into 
the diffuser and that the flow in the core is one-dimensional and isen- 
tropic. Although this model can be taken only as a rough approximation 
of the physical case, it was hoped that its results would be useful at 
least for the prediction of trends* A brief description of this method 
follows . 


The method uses a momentum integral approach to the calculation 
of the boundary layer characteristics* The velocity at the wall is allowed 
to be finite to include the effects of suction* Definitions of momentum 
and displacement thicknesses which account for transverse curvature are 
used. The solution for the momentum integral equation is 

X 

®co “ J^w(x) G(x) dx, (14) 

o 
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where 


G(x) = exp 


r(5*/e) 


u 


ad du ^ 


dx a p dx 

*^00 00 


(p a) 

' >^00 00 


dx 


(15) 


and 


C- V T 

r w 00 

W(x) = ~ sec ^ ~ 

00 W 


(16) 


Flat-plate values for 5*/ 0 and Cf were used in the above equa- 
tions, These were obtained from an exact solution developed by Iglisch 
(1949) for incompressible flow. In Iglisch* s work, these quantities are 
functions of the parameter 


V 


w 


u 

00 


n/rc 


inc • 


(17) 


But since 


Re. = (Re ) = (T (Re^^ ); 

me wall 00 w comp 


(18) 


the results of Iglisch can be used directly as functions of 


V 

w 

u 

00 



com] 


py 


(19) 


to obtain Cf and 


(6*/e) 


2 d inc* 


The value of (S*/e) 2 d inc 
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in the following equation (Sivells and Payne, 1959) to obtain 


(s*/e) 


2d 


(6*/e) 


T T 

w ^ aw 

2 d inc T T 

00 00 


1 . 


( 20 ) 


The parameter v^/u was computed by assuming that the flow through the 
porous wall is choked with stagnation conditions equal to the wall static 
pressure and temperature. The hole flow coefficients were taken from 
experimental data for thin orifices over a Reynolds number range from 
continuum to near free molecule. 

Equation (14) must be solved numerically. A program for an 
IBM 7090 computer was developed and was available for the present study. 
No account was taken of pressure gradient effects on Cf and 5*/ 0. 


B. Modification of Boundary Laver Equations to Include Pressure 
Gradient Effects 


During the present study, an effort was made to include the 
effects of pressure gradient in the IBM 7090 computer program, because 
of the possible importance of the adverse pressure gradient in a diffuser. 

No solutions were found in the literature for the compressible 
laminar boundary layer with heat transfer, suction, and pressure gradient. 
Thus, to estimate the pressure gradient effect, it is necessary to resort 
to solutions of the solid wall boundary layer with pressure gradient and 
to apply these as perturbation type corrections to the Cf and 6*/ 0 values 
for the suction boundary layer. 

Solutions to the laminar boundary layer equaLions that Include 
the effects of heat transfer, compressibility and pressure gradient have 
been quite limited in number. The method of Cohen and Reshotko (NACA- 
TR-1294, 1956), as distinct from other methods, does not require the 
solution of one or more ordinary differential equations and seemed to 
be the most suitable for inclusion into the method of Bottorff and Rogers. 
Cohen and Reshotko, after applying Stewartson’s transformation to the 
boundary layer equations, use Thwaite*s concept (developed by him for 
incompressible flow) of relating the wall shear, its normal derivative 
at the wall, and the form factor to one another without specifying a 
type of velocity profile. Nondimens ional forms of these quantities 
were defined and were evaluated by examining exact solutions for the 
laminar boundary layer, in this case those of Cohen and Reshotko (NACA 
TR-1293, 1956). 
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The Cf equation developed by Cohen and Reshotko is the following: 


( 21 ) 


Converting this to the Cf used in Bottorff and Rogers, we have 


( 22 ) 


Cohen and Reshotko have obtained the shear parameter, as a 
function of the correlation number, np, which must be known if this 
method is to be used. If there is no pressure gradient, np = 0. 

Since the Cohen method was not developed for a suction boundary 
layer, application of its results to a porous nozzle-diffuser is limited 
to a "perturbation” type of correction to the Cf developed in Bottorff 
and Rogers for suction boundary layers. Thus, a Z^f can be defined as 





(23) 


where is the pressure gradient correction to the shear parameter, 
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In the case of 5*/0, 

equation: 


Cohen and Reshotko give the following 


H = ^ = + 

9 tr 


(H + 1), 

2 00 ' tr 


( 24 ) 


where Hj-r is a transformed form factor and is, like Cj, a function of 
the correlation number, np. For no pressure gradient. 


«tr = «inc T 


w 


aw 


If a pressure gradient exists, can be written 


«tr = “inc 

aw 


(25) 


where Z^tr correction due to pressure gradient from Cohen and 

Reshotko. If Htr is used instead of Hj_nc the method of 

Bottorff and Rogers, the resulting equation for (6*/0)2d with pressure 
gradients is 


(s'Ve) 


2d 



+ 


w 


tr 




(26) 


C. Considerations Regarding Separation 

In modifying the computer program to include diffuser calcula- 
tions, the range of the pressure gradient correlation parameter, np, had 
to be arbitrarily limited to a value below that which would cause separa- 
tion for a no-suction boundary layer. The possible error introduced by 
doing this may be large, as indicated by Figure 3, where a typical cal- 
culated variation of np along the diffuser is plotted. It is clear that, 
if the boundary layer actually separates at the point indicated for no 
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suction, diffuser performance will be poor. Although the plot indicates 
that a no-suction diffuser would separate at the entrance, the np at this 
point is not considered valid because it occurs at the nozzle-diffuser 
boundary. If the boundary layer separates at even the second indicated 
point, however, the contraction of the diffuser would be quite limited 
and the pressure recovery would be essentially that associated with a 
test section normal shock. It has been found experimentally, however, 
that this separation apparently does not occur, or if it does occur, it 
fails to seriously disrupt the core flow. This statement is based on 
the fact that pressure recoveries of up to ten times normal shock were 
achieved in the tests and also upon the fact that measured diffuser static 
pressures exhibited a smooth increase along the length of the diffuser. 

The separation point indicated on the plot was for a no-suction 
boundary layer. It is well known that boundary layer suction delays 
separation, so that a higher value of the separation np can be expected 
for a suction boundary layer. Investigations have been made to attempt 
to determine the magnitude of the increase in the pressure gradient 
parameter that will be caused by the suction. No references were found 
relating the effect of suction to the separation pressure gradient param- 
eter in a compressible boundary layer. In view of this lack of informa- 
tion on compressible boundary layers, it was decided to attempt to 
determine the order of magnitude change in pressure gradient parameter 
due to suction by using the results of incompressible analysis. 

Spalding and Evans have prepared a series of reports that com- 
pile the available exact solutions of the incompressible boundary layer 
with an arbitrary pressure gradient and suction. In their work, the 
suction boundary layer is characterized by parameters such as the momen- 
tum thickness, the kinematic viscosity, the local axial velocity gradient, 
and the suction flow rate. The velocity gradient is assumed to be of the 
form du/dx = cu^. The sign of the constant, c, will be positive or nega- 
tive, depending upon whether the flow is accelerating or decelerating. 

For all of the data presented, the exponent, n, has been limited to 
values less than 2 for positive velocity gradients and to values greater 
than 2 for negative velocity gradients (diffusers). The limitations 
have no physical significance, but are made only for mathematical simplic- 
ity. These restrictions are important in the case of hypersonic diffusers, 
however, since for the velocities of interest in the present study (3000 
to 4000 fps) calculations have been carried out only for extreme velocity 
gradients. It is thus apparent that the incompressible flow calculations 
are not useful in attempting to determine the magnitude of the effect of 
suction on the separation point. 
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D, Results of Calculations With and Without Pressure Gradient 


Figure 4 is a comparison of the computer results with and with- 
out the previously described pressure gradient correction. Mach number 
and m^/m^ have been plotted versus the distance along the nozzle-diffuser 
centerline for the 8-degree half-angle, 15 percent porosity diffuser at 
stagnation conditions typical of those used in the tests. The pressure 
gradient correction gives a lower Mach number at all locations; this 
implies that the increased skin friction coefficient more than offsets 
the favorable shape parameters and results in an increase in boundary 
layer thickness. However, the Mach number difference seems to disappear 
near the end of the diffuser. Mass flow ratios, however, are about 5 to 
15 percent higher (suction flow 5 to 15 percent lower) in the no-pressure 
gradient case, a result of the lower tunnel static pressures for this 
case. 


To further assess the validity of the Cohen-Reshotko method for 
calculating laminar boundary layer characteristics in a pressure gradient, 
the nozzle boundary layer was calculated for the pressure gradient and 
no-pressure gradient cases for an existing M = 6 nozzle, which has been 
thoroughly investigated experimentally. The comparison between the 
resulting test section Mach numbers is shown in Figure 5. Again, the 
pressure gradient case shows a decreased Mach number and is in poorer 
agreement with experiment than the no-pressure gradient calculations. 

It appears that the Cohen-Reshotko procedure may over-estimate 
the effect of the pressure gradient on the skin friction coefficient. 

This is also implied in a report by Carden who compared experimentally 
measured heat transfer coefficients in a laminar flow nozzle with the 
results of calculations using the Cohen-Reshotko procedure. The heat 
transfer coefficients calculated by the Cohen-Reshotko procedure greatly 
exceed experimentally determined ones. In view of this uncertainty of 
tbft adequacy using the Cohen-Reshotko procedure to calculate the effect 
of pressure gradient on noz^le-diffu?pr boundary layer characteristics 
and also considering that the effect of the pressure gradient appears 
small when used in the calculation of the nozzle-diffuser boundary layer, 
for the purposes of this study the effects of pressure gradient can be 
neglected. 


IV. RESULTS OF THEORETICAL CALCULATIONS 


A. Mass Flow 

Calculations were made for a series of nozzle-diffuser combina- 
tions. The variables included the Mach number, the Reynolds number, the 
contraction angle, w, and the porosity, cp. While the important parameters 
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in determining the performance of a diffuser are the pressure recovery 
and mass flow relationships, by the use of equations (9) and (10), it is 
possible to determine the limits on pressure recovery if the variation of 
mass flow and temperature ratio with area ratio is known. Figures 6 
through 8 show the typical variation of mass flow ratio with diffuser 
area ratio. 

The porous nozzle used in these calculations had a 26-degree 
included angle and a 12-inch exit diameter. The nozzle throat was 
varied between Mach number 10 and Mach number 6, The lowest values of 
A^/As are those corresponding to the test section. 

Figure 6 shows the variation of mass flow with area ratio for 
different values of cp and oo. This plot shows that increasing either the 
porosity or the length increases the flow through the pores. This is 
consistent with what would be expected intuitively, i.e., that the flow 
out the diffuser wall would be approximately proportional to the total 
open area of the diffuser. The total flow through the diffuser pores is 
not directly proportional to the area of the diffuser pores because 
increasing diffuser pore flow results in decreasing static pressures in 
the diffuser. Since the local mass flow through the diffuser pore is 
almost directly proportional to the local static pressure, the decrease 
in static pressure results in a decrease in mass flow per unit open area. 
This assumed relationship between local mass flow and local static pres- 
sure keeps the flow from being directly proportional to the total open 
area of the diffuser. 

Figures 7 and 8 show the effect of Reynolds number on the dif- 
fuser performance. Increasing Reynolds number results in a relative 
reduction in mass flow through the pores. This is due to the increase 
in Mach number that is associated with an increase in Reynolds number. 

The Mach number increases because of the reduction in boundary layer 
height with increasing Reynolds number. The increased Mach number 
results in a relatively lower static pressure at the pores, and since 
the pore flow is directly proportional to the static pressure, this 
results in a reduction in pore flow. As the Reynolds number increases, 
the boundary layer becomes thinner, so that the variations in the height 
of the boundary layer have a smaller effect on the Mach number and pres- 
sure in the nozzle and diffuser. 

B. Heat Transfer 


To use Figures 6 through 8 in calculating the pressure recovery, 
it is necessary to determine the variation of Calculations 

were made using the procedure outlined on pages 7 and 8. It was found 
that the heat transferred between the nozzle throat and the diffuser 
throat was relatively constant and equal to about 15 percent of the 
energy entering the nozzle throat. Thus, the temperature ratio (T^p^/T^j-^^) 
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is about .85, Since only a small fraction of the energy is transferred 
from the incoming gas to the nozzle walls, it is apparent that the tem- 
perature ratio is not an important parameter in calculating the theoret- 
ical limits on the pressure recovery. 

C. Cryopumped Diffusers 

In a low density wind tunnel with porous diffuser, the flow 
through the pores generally represents a very large volume flow rate 
because of the low pressures (typically of the order of 1 micron) which 
must be maintained outside the nozzle and diffuser walls. Thus, a cryo- 
pump, which can be arranged so that it entirely surrounds the nozzle and 
diffuser, is an especially attractive means of pumping the suction flow 
because of the very high pumping speeds which can be achieved. 

Since the cryopump is, however, a mass flow limited pump, the 
problem in such an installation is to design a diffuser which will mini- 
mize the volume flow rate at the diffuser throat while maintaining the 
pore mass flow below the capacity of the cryopump. To illustrate the 
trade-off involved, calculations were made for the 12-degree half-angle, 
15 percent porosity diffuser with the 12-inch diameter Mach 10 nozzle. 
The test section unit Reynolds number was 340/ in (P-p = 1.96 psia, 

TTq = 1320°R), Figure 9 presents the maximum possible pressure recovery 
(equation 10) and the diffuser throat mass flows as a function of nozzle 
to diffuser throat area ratio. As the diffuser area ratio is increased, 
the assumed normal shock occurs at a lower Mach number. The lowered 
Mach number more than offsets the decreased mass flow ratio to give a 
steadily increasing pressure recovery as the area ratio is increased. 
Figure 10 presents the corresponding pumping speed requirements for the 
diffuser throat flow along with the cryopump capacity required to pump 
the pore flow. The diffuser throat volume flow can be diminished to 
very small values, but the cryopump capacity curve rises steeply as this 
is done . 


It is apparent then that a trade-off musL be made fn balance the 
two requirements. The most economical contraction ratio cannot be derived 
in general terms since it depends on a great many factors (Reynolds number 
range, tunnel size, availability of various pumping means, etc.) which 
must be individually considered for each installation. 

D. Starting Calculations 

The porous diffuser can start with greater amounts of contrac- 
tion, than a solid wall diffuser since the openings in the diffuser wall 
downstream of the shock system act as additional diffuser throat area 
during the starting process. Since these pores are not carefully shaped 
nozzles, but are sharp-edged orifices, instead of flowing full of gas at 
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a sonic condition, they flow at a lower rate. While, if the pressure 
ratio is high enough, the flow rate can be estimated as a function of 
Reynolds number, because of other uncertainties, it is adequate to 
assume that the pores are choked. 

The temperature of the gas passing through the pores is unknown 
since it is difficult to estimate the heat transfer. It was found that 
the heat transfer in the diffuser is low if the diffuser is supersonic 

^ -85); however, if the gas is subsonic and at a higher pressure 
level, the heat transfer will be increased. In view of this uncertainty, 
a conservative assumption has been made; i.e., the gas is all at the 
stagnation temperature. 

The other major uncertainty is the pressure recovery of the gas 
behind the normal shock system. At high Reynolds number, it has been 
found experimentally that this shock system has a lower loss than a single 
normal shock. This is evidenced by the ability of diffusers to start with 
greater contraction than theory predicts. For the present analysis it has 
been assumed that these errors are self-compensating; i.e., the reduction 
in area is just balanced by the higher pressure recovery and lower tem- 
perature. Now the starting criterion will be that the quantity of gas 
that can flow through the open area downstream of the shock wave must 
equal or exceed the quantity of gas entering the shock wave. The gas 
flowing through the pores is assumed to be choked at the stagnation tem- 
perature and pressure associated with a normal shock. The starting 
characteristics of several nozzle-diffuser configurations were investi- 
gated using this approach. It was found that the crucial phase of the 
starting process occurs when the normal shock is positioned at the test 
section. This is the same result found for a nonporous diffuser. 


V. EXPERIMENTAL PROGRAM 


A. Facility Description 

The experiments were conducted in the Hyperaltitude Facility of 
the Environmental Division of the U. S. Naval Missile Center, Point Mugu, 
California. A complete description of this facility is given in Bottorff 
(1964), and only a brief summary will be presented in this section. The 
basic facility consists of a 10-foot diameter, 20-foot long vacuum chamber 
that has a combination of pumping systems. The primary pumping unit for 
the present test was a 20°K cryopump that is cooled by a 350 watt gaseous 
helium refrigerator. The condenser can be isolated from the main chamber 
by a large 5-foot diameter valve which allows access to the models and 
tunnel without bringing the cryopump up to ambient temperature. 
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The present investigation used the porous wall nozzle described 
in Bottorff and Rogers (1963), This nozzle is basically a Mach 10 nozzle 
with 12- inch exit diameter which has liquid-nitrogen-cooled porous walls. 

B. Experimental Arrangement 

The experimental arrangement used in this test is illustrated 
in Figure 11, The nitrogen gas which was obtained by vaporizing liquid 
nitrogen was metered into the stagnation chamber to maintain a selected 
stagnation pressure. After passing through the nozzle throat, part of 
the flow was removed through the liquid-nitrogen-cooled porous walls of 
the nozzle and the diffuser. This flow passed through the 5-foot diameter 
valve into the cryopump. The remaining flow passed through the diffuser 
throat into the pumping system where it could be directed to either the 
mechanical pump or the cryopump. The pressure at the end of the diffuser 
was varied by manipulating the throttle valves (1 and 2), 

It was originally planned to use the section M (Figure 11) as a 
metering run; however, this did not prove practical because of erratic 
pressure drops in the metering section. This erratic performance resulted 
from the persistence of the core flow when the diffuser throat was super- 
sonic, This problem was alleviated by installing the stilling chamber S 
(Figure 11) and using the thin-walled orifice as a metering system. To 
use this metering system, it was necessary that all of the flow pass 
through the stilling chamber and into the cryopump so valve 2 was closed 
at all times, 

C. Instrumentation 


Since there are extremely wide variations of pressure level in 
the hypersonic wind tunnel, it is necessary to use a variety of pressure 
gauges to monitor the flow conditions, 

me s L.agna tion chamber pressure was measured using a Bourdon 
gauge which was limited by reading accuracy to 3 Lo 5 percent- accuracy. 
The static pressures in the converging section of the supersonic diffuser 
were measured using thermocouple gauges. The first two instruments had a 
usable range of 0-100 microns, while the third gauge had a range of 
5-1000 microns. By the use of suitable valving, the test section noirmal 
shock pressure, the static pressure in the diffuser throat, and the dif- 
fuser recovery pressure were measured using one Alpha tron gauge. When it 
was necessary to measure the static pressure near the diffuser recovery 
tube, the same Alphatron was used. The pressure in the metering system 
was measured using a capacitance type mechanical diaphragm gauge. All 
gauges were calibrated on a device that metered known increments of gas 
into a container of a fixed and known volume. 
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D. Diffuser Configuration 


Three diffuser contraction sections were tested. Two contrac- 
tions had an 8-degree half-angle and the third had a 12-degree half-angle. 
One 8-degree diffuser and the 12-degree diffuser had a porosity of 15 per- 
cent; i.e., there was 15 percent open area. The second 8-degree diffuser 
had a porosity of 30 percent. The pores consisted of 1/4- inch diameter 
holes drilled in the 1/16-inch diffuser wall. The diffusers were con- 
structed so that the diameter of the throat section could be varied 
between 2 inches and 31/2 inches. The L/D of the throat was essentially 
constant at a value of 4.8. The entire diffuser assembly could be moved 
relative to the nozzle so that the effect of varying the free-jet length 
could be investigated. The supersonic contraction section and the con- 
stant area throat were cooled to liquid-nitrogen temperature. 

The configurations are identified by a code consisting of dif- 
fuser half-angle, porosity, throat diameter and free-jet length. Thus 
8®-157o-2”-6" refers to the 8-degree half-angle, 15 percent porosity 
diffuser with a 2-inch diameter throat and a free-jet length of 6-inches. 

A traverse mechanism was installed so that various models could 
be inserted in the stream. Since these were expected to be high per- 
formance diffusers, it was anticipated that the disturbances caused by 
the models would seriously affect the diffuser performance. The models 
included a 1-inch diameter sphere, a 30-degree included angle cone with 
a 1-inch base diameter, and a 1/4- inch stagnation pressure probe. 

E. Testing Procedure 

During typical testing, the tunnel conditions were established 
by raising the stagnation chamber pressure to the desired level with the 
throttle valve wide open. After the pressures had stabilized, the read- 
ings were taken, and the throttle was closed to raise the pressure at 
the end of the diffuser section. Typically the valve was closed until 
the metering system showed a decrease in mass flow, and then the valve 
was opened until the entire mass flow was once more passing through the 
metering section. This point would roughly correspond to the critical 
point of the diffuser. Additional data points were then taken with the 
valve closed beyond this point. This procedure generated plots of pres- 
sure recovery versus mass flow that are comparable to those obtained 
during the testing of supersonic inlets. 

It was found that some configurations would not operate properly 
because of the losses in the 4- inch diameter metering piping. Since the 
mass flow data were the most important, the metering system was shortened 
and the valves were removed in an attempt to reduce the losses through 
the metering system. For these configurations, only a single mass flow 
point was obtained for each Reynolds number. 
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F. Data Reduction 


The pressure data were reduced to coefficient form by ratioing 
the pressures to the stagnation chamber pressure. 

The mass flow through the metering system was obtained by assum- 
ing that the flow characteristics of the orifice were the same as those 
reported by Leipmann (1960). In that investigation, an orifice with an 
L/D of 1/40 was tested from the continuum range through the free-molecular 
flow range. Liepmann found that if the Knudsen number of the orifice was 
below .1, the mass flow was essentially independent of Knudsen number and 
equal to 85 percent of the flow that would pass through a sonic throat of 
the same area. The present experiments were in the same Knudsen number 
range; however, it was not possible to maintain the same pressure ratio 
across the orifice. Liepmann maintained a pressure ratio of 1,000. How- 
ever, he points out in the theoretical development that in continuum flow 
a pressure ratio of approximately 26 is sufficient to prevent the down- 
stream pressure from influencing the flow through the orifice if the gas 
has a value of 7 = 1.4. This condition was met for most of the configura- 
tions tested in the present investigation. 

The .relationship between the flow through the nozzle and the 
flow through the diffuser throat is given by 



o 


In the present investigation, it was found that, with room temperature 
gas entering the nozzle, the temperature of the gas leaving the orifice 
was also room Leiuperature. The orifice diameter was 4 inches and the 
nozzle throat diameter was 9/32 inch. Equation (27) becomes 



(28) 


This is the expression that was used in calculating the diffuser throat 
mass flow ratio from the measured pressures. 
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VI, EXPERIMENTAL RESULTS AND COMPARISON WITH THEORY 


A. Pressure Recovery Data 

Figure 12 is a plot of measured pressure ratio versus measured 
mass flow ratio for the diffuser configuration 8°-15%-2”-0. This fig- 
ure can be used to demonstrate some of these salient features of the 
experimental results. 

The flagged symbols correspond to the pressure recovery meas- 
ured just downstream of that location. The numbers on the points corre- 
spond to the sequence of closing the mass flow valve. While the static 
pressure measured at the aft end of the diffuser follows the theoretical 
trend insofar as the pressure rises as the valve is closed, the stagna- 
tion pressure often does not. This is because the single stagnation pres- 
sure tube is generally not representative of the flow conditions at that 
station* Since the probe is mounted in the center of the tube, it is 
unduly influenced by the high stagnation pressure core which persists 
into this region when the throttle valve is relatively wide open. As 
the valve is closed, the shock system is moved forward and promotes 
better mixing; thus, the stagnation pressure probe becomes more repre- 
sentative of the flow at that location. 

As the shock system is moved forward the diffuser mass flow 
begins to decrease. This is probably the result of separation and 
reversed flow on the walls of the diffuser throat. This separation 
feeds far enough forward to influence the static pressure measurements 
in the converging section of the supersonic diffuser. 

Because the disturbance feeds forward such a long distance, 
the peak pressure recovery, is not obtained until a considerable part of 
the mass flow is being spilled at the low Reynolds numbers. While this 
could probably be corrected by using a longer constant area throat, the 
main benefit would be the increased mass flow at the critical point, 
since the increase in pressure recovery would be modest. This is clear 
from the comparison of the experimental pressure recovery with the theo- 
retical maximum possible pressure recovery. The maximum theoretical 
pressure recovery is obtained by the intersection of vertical lines 
through the maximum measured mass flows and the line labeled equation (10). 
This corresponds to isentropic flow between the nozzle throat and the dif- 
fuser throat, and therefore represents the maximum possible pressure 
recovery. The maximum experimental pressure recoveries are within 10 to 
15 percent of the theoretical values. 
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The agreement between the values of pressure recovery suggests 
that the relatively short constant area diffuser throats are adequate if 
some spillage is allowable. If the diffuser walls were solid so that 
mass flow could not be spilled without causing the diffuser to become 
unstarted, the pressure recovery would be significantly reduced. Under 
these conditions, the maximum pressure recovery would be limited to the 
values corresponding to the point at which the diffuser throat mass flow 
starts to decrease. 

Figure 13 illustrates the change in diffuser effectiveness 
caused by opening the free-jet to 6 inches and by installing a model. 

To avoid the erratic pressures measured by the stagnation pressure probe, 
the pressure ratios were obtained using the static pressure at Station 3. 
For most configurations, this pressure is within 10 percent of the stag- 
nation pressure at the peak pressure recovery. 

The 6-inch free-jet causes a slight reduction of the mass flow 
and pressure recovery. The model has a more deleterious effect on the 
performance, causing significant reductions in mass flow and pressure 
recovery. 


Figures 14 and 15 compare the same performance parameters for 
configurations 8°-30%-2 and 12°-15%-2,5, respectively. These configura- 
tions show similar characteristics to the ones previously discussed. 

It was not always possible to obtain continuous data as the 
mass flow was reduced by closing valve 1, As the mass flow was reduced, 
a point was often reached where the shock wave would move abruptly from 
near the diffuser throat to a point far upstream of the diffuser throat. 
This would result in negligible values of pressure recovery and mass 
flow. It appears that, under these circumstances, the suction was insuf- 
ficient to stabilize the shock wave system at that point in the diffuser. 

Similarly, it was not always possiMe to obtain data at the 
higher Reynolds numbers. For example, the configuration 12"-i5/>3**-G** 
would not remain started above Reynolds number/inch = 800. As the stag- 
nation pressure was raised above this point, the shock system would 
abruptly move upstream from the diffuser throat. Since this configura- 
tion had a pressure recovery that was comparable to the pressure required 
to overcome the losses in the metering system, it is possible that the 
metering system losses were acting as a partially closed valve and forc- 
ing the shock system forward of the diffuser throat. In this case, the 
same mechanism would be involved in the abrupt movement of the shock 
wave. In one case, the crucial condition is achieved by varying the 
mass flow throttle, while in the other case, it is achieved by raising 
the stagnation pressure. 
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From a comparison of the experimentally determined variation 
of pressure recovery with mass flow and the theory given by equations 
(9) and (10), it is apparent that the theory serves as adequate limits 
to diffuser performance. 

B. Mass Flow Data 


Since the data presented so far have confirmed the validity of 
the limits set by equations (9) and (10), the remainder of the discussion 
will be concerned with mass flow ratios. Figures 16 through 19 illus- 
trate the variation of mass flow ratio with Reynolds number for the 
various configurations. The theoretical test section mass flow is 
included to indicate the magnitude of the flow through the nozzle pores. 

Figures 16, 17, and 18 compare the mass flow characteristics 
of a series of configurations which only differ by the diameter of the 
constant area throat. The 2-inch diameter throat and the 3-inch diam- 
eter throat both follow the theoretical trend, while the 3 1/2-inch 
diameter throat has a lower mass flow than expected. This reduced 
mass flow was caused by excessive losses in the mass flow metering 
system which prevented the 3 1/2-inch configuration from starting prop- 
erly, This is illustrated by Figure 20, which shows the measured 
static pressures in the converging section of the diffuser for the 
three configurations. The 3 1/2-inch configuration shows a distinct 
rise in the static pressure before the throat. This pressure rise 
increases the flow through the pores and results in a lowered mass flow 
at the diffuser throat. 

Figures 16 and 19 are the same configuration except for an 
increase in porosity. The increased porosity had a small effect on 
the flow rate in the diffuser. This is in contrast to the theory, 
which predicted a significant increase in flow rate through the pores. 

It was found that the configuration with increased porosity had a lower 
static pressure distribution along the converging section of the dif- 
fuser. This lowered static pressure would yield a reduced mass flow 
rate compared to the 15 percent configuration. It will be shown in a 
later section that, while the level of the static pressure does not 
appear to be a satisfactory indication of the flow through the pores, 
it appears that the variation of the static pressure does coincide 
with the variation of mass flow. This is borne out by the changes 
which occurred when a model was introduced into the stream. The 
introduction of the model caused very high flow rates through the pores, 
and a significant increase in static pressure was measured along the 
converging section of the diffuser. 

In general, the theoretical and experimental flow rates are in 
good agreement for the configurations with no free-jet. The six-inch 
length of free-jet reduced the mass flow by 5 to 10 percent for most 
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configurations. The introduction of the model caused a significant 
reduction in mass flow at all Reynolds numbers, but the influence was 
most pronounced at the low Reynolds numbers. 

The agreement between the theoretical and the experimental 
values of mass flow must be regarded as somewhat fortuitous, since the 
measured static pressures in the contraction section of the diffuser 
were significantly above the theoretical values in all cases. Figure 
20 is typical of this comparison. In the theoretical analysis, the out- 
flow through the pores was assumed to be nearly proportional to the local 
static pressure. It would be expected that the increased static pressure 
would result in a much higher pore mass flow rate than that predicted by 
the theory, which used lower static pressures. This would reduce the 
mass flow rate through the diffuser throat compared to the theory. For 
example, if a calculation is made of the pore flow using the experimen- 
tally measured pressures (Figure 20) instead of the theoretical values, 
the mass flow through the pores would be more than doubled. This did 
not occur experimentally, and no explanation has been found for this 
behavior , 


A limited investigation was made of the effect of varying the 
chamber pressure. It was found that some configurations were unusually 
sensitive to small chamber pressure variations. For example, during 
typical operation, the chamber pressure varied between 1/2 and 2 microns, 
depending upon the flow rate. The pressure was limited by the conduc- 
tance of the 5-foot diameter valve which leads to the cryopump. It was 
found that raising the chamber pressure by 1 micron could cause a signi- 
ficant change in mass flow and pressure recovery. This is illustrated 
in Figure 21 which shows the variation in diffuser throat mass flow with 
chamber pressure for a configuration with and without a model installed. 
Since the pressures measured on the inside of the diffuser were of the 
order of 10 microns or greater, it is difficult to understand how a 
chamber pressure variation of 10 percent of this value would cause such 
a large change in mas« flow. It seems possible that the phenomenon that 
is responsible for this variation of pore flow with chamHpr pressure may 
also be responsible for the reduced values of pore flow when compared to 
the theoretical values associated with the high static pressure levels 
measured in the diffuser. 


VII. COMPARISON WITH PREVIOUS EXPERIMENTS 


Rogers (1962) presented the results of an experimental investigation 
of a porous cooled diffuser. In these experiments an attempt had been 
made to fabricate a variable porosity diffuser. This was done using two 
porous concentric cones. Because of the difficulty in maintaining align- 
ment between the pores, the level of porosity was not accurately established. 
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A further difficulty was encountered because only the outer cone was 
cooled with liquid nitrogen. The cooling on the inner cone depended 
upon conduction from the outer cone. In the present analysis, it has 
been assumed that the wall was cooled to liquid nitrogen temperatures. 

The effective porosity was determined by making calculations for a 
series of porosities and selecting the one that agreed with the experi- 
ments at a specific Reynolds number. It was found that, whereas the 
porosity of the individual cones was 30 percent, the effective porosity 
was 7,5 percent. The comparison between the theoretical and experimental 
variation of mass flow with Reynolds number is shown in Figure 22. The 
porosity had been selected to agree with the experiment at Re/ in = 1000. 

It is seen that there is general agreement at all Reynolds numbers except 
the highest. This highest Reynolds number point is somewhat questionable 
since the flow leaving the nozzle exit was nonuniform and had strong 
compressions on the outer edge. While it would be expected that the 
strong compressions would raise the static pressure and consequently the 
flow out through the pores of the diffuser, this does not appear to have 
been the case. The mass flow passing through the diffuser throat was 
greater than that predicted by the theory, which indicates a lower mass 
flow passing through the pores. This is another example of the difficulty 
of calculating the flow through the pores. 

Figure 23 shows a comparison between the theoretical and experi- 
mental values of pressure recovery. In this case, equations (9) and (10) 
have been used to calculate the limits on the pressure recovery, but the 
experimentally measured mass flow has been used in place of the theoretical 
value. The results show that the experimental pressure recoveries are 
slightly below the minimum theoretical values. This is the result of 
losses occurring downstream of the throat of the diffuser. 

This diffuser configuration was also adversely affected by the 
presence of the model in the test section. When a cylinder was inserted 
into the stream so that it spanned the tunnel, the mass flow was decreased 
by approximately one-half. The cylinder had a diameter equal to about 
6 percent of the test section diameter. 


CONCLUSIONS 


The results of the series of investigations on porous wall low den- 
sity wind tunnel diffusers have shown good agreement between theory and 
experiments. It was predicted theoretically and confirmed experimentally 
that the limits on the pressure recovery at the diffuser throat can be 
given by the following simple equation: 




n/TtTtT 


(29) 
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In this expression, the constant, C, is unity as a lower limit and 1.65 
as an upper limit for hypersonic flow with 7 = 1.4. The lower limit 
corresponds to choking at the diffuser throat, while the upper limit 
corresponds to isentropic flow between the nozzle throat and the dif- 
fuser throat followed by a normal shock. From an analysis of several 
nozzle-diffuser combinations, it was concluded that the heat transferred 
to the walls was small and the term would not vary significantly 

from .85. 

A comparison of the results of a computer program and the experi- 
mental investigation showed that the theory accurately predicted the 
variation of diffuser throat mass flow with diffuser throat area ratio. 

It was concluded, however, that this agreement was somewhat fortuitous 
since the level of static pressure in the diffuser was much higher than 
the theoretical value. This higher static pressure should have resulted 
in an increased mass flow out through the porous walls of the diffuser 
and a resultant decrease in mass flow through the diffuser throat. Since 
this was not found to be the case, it was concluded that the model used 
in calculating the flow through the diffuser pores was incorrect. It is 
suggested that this is an area that warrants further study. 

It was found that introducing a model into the stream caused a 
significant reduction in diffuser throat mass flow and pressure recov- 
ery. The large decrease in diffuser throat mass flow implies an increase 
in mass flow through the porous walls. The diffuser static pressures 
increased when the model was installed. This would lead to an increase 
in flow through the porous walls if the flow is proportional to the static 
pressure. 

The theoretical investigation indicated that it was possible to 
obtain pressure recoveries approaching unity for the flow remaining in 
the diffuser throat. These high pressure recoveries could be obtained 
only at the expense of having most of the flow pass through the porous 
walls ot the uiffucer. Since any practical design must consider the 
quantity of flow passing through the ditfuser porcc as 1 as the flow 
passing through the throat, it is necessary to make a trade-off between 
diffuser throat pressure recovery and pore mass flow. 
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Figure 1. Wind Tunnel Schematic 
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Figure 5. Comparison Between Theoretical 
and Experimental Mach Numbers, 





Figure 6. Variation of Diffuser Mass Flow 


















CRYOPUMP 



Figure 11. Experimental Arrangement 
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Figure 16 . Variation of Mass Flow with Reynolds No 
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Figure 17 . Variation of Mass Flow with Reynolds No 
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Figure 18. Variation of Mass Flow with Reynolds No 
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Figure 25 - Comparison of Theoretical 
and Experimental Pressure 
Recovery 
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